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Abstract 

For OFDMA systems, we find a rough but easily computed upper 
bound for the probability of loosing communications by insufficient num- 
ber of sub-channels on downlink. We consider as random the positions of 
receiving users in the system as well as the number of sub-channels dedi- 
cated to each one. We use recent results of the theory of point processes 
which reduce our calculations to that of the first and second moments of 
the total required number of sub-carriers. 

1 Introduction 

The demand for high data rate wireless appHcations with restrictions in the RF 
signal bandwidth requires bandwidth efficient air interface schemes. It is known 
that OFDM yields a relatively simple solution to these problems [Tj. Based on 
the OFDM system, OFDMA can achieve a larger capacity. Furthermore, this 
latter system is more flexible, since it can be easily scaled to flt in a certain piece 
of spectrum simply by changing the number of used subcarriers However, 
as any wireless systems, OFDM and OFDMA have physical limitations which 
cause loss of communications. This loss can be caused by insufficient power or by 
low signal-to-interference ratio, for instance. In this paper we are interested in 
the calculation of an upper bound of the probability of loosing a communication 
due to an insufficient number of sub-channels in the downlink. 

We say that the system is overloaded when all non-used sub-channels are 
not enough to warrant a minimum data rate for an incoming demand. We 
consider a system with Nq sub-carriers and Ni is the number of sub-carriers 
used by the i-th user in the cell. As it is usually done, we substitute the 
flnite number of subcarriers by inflnity and substitute the loss probability by 
Pioss — P{J2i^i > ^o)- It is well known that this consideration gives us an 
upper bound for the actual loss probability. 

A user i requires a capacity C, depending on the service he uses. Considering 
a system with just one kind of service, all users require the same capacity Cq. 
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Even so, the number of subcarriers for each user varies according to the channel 
conditions. These conditions can be summarized into two kinds of gains, one 
depending only on the position of the user z, the path loss Gp/. , and a gain 
Gi, which classically may include the shadowing and the Rayleigh fading. We 
choose the simplest model to represent the path loss: 

i 

where K and 7 are constants and is the distance between the user i and the 
antenna. 

Shannon's maximum achievable capacity implies that: 



M/log,(l + ^^^) 

where W is the bandwidth of each sub-carrier, Pt is the mean transmitted power 
by sub-channel and / is the power of the additive Gaussian white noise by sub- 
channel. 

We consider that the number as well as the position of users in the cells 
are random. After some natural assumptions done in the following section, we 
conclude that the configuration of users in the cell is a Poisson point process 
(see section [2]). 

After a summary on Poisson point process, we consider three different cases 
to calculate an upper bound for the loss probability. First we consider the 
simplest case with deterministic gain. In Section lU we consider a non-selective 
frequency gain, the shadowing. In section O we consider a general case from 
which all other cases could be derived but for which no closed form formula 
exists. 



2 Poisson point processes 

For details on point processes, we refer to [H, 0, S S|- ^ configuration 77 in R'^ 
is a set {x„, n > 1} where for each n > 1, x„ G R'^, Xn 7^ Xm for n ^ m and 
each compact subset of R'^ contains only a finite subset of 77. We denote by V-Rk 
the set of configurations in R'^. Equipped with the vague topology of discrete 
measures, V-^k is a complete, separable metric space. A point process <& is a 
random variable with values in Fp^fc, i.e., = {Ar„(tj), rt > 1} S Fp^fc. For 

A C R*^, we denote by the random variable which counts the number of 
atoms of ^(lS) in A: 

= ^lx„MeA GNU 
Poisson point processes are particular instances of point processes such that: 
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Definition 1. Let A be a a finite measure on R'^. A point process $ is a 
Poisson process of intensity A whenever the two following properties hold. 

1 - For any compact subset A e R'^, $a is a Poisson random variable of 

parameter K(A), i.e., 

P(<I>^ = fc) = e-^(^)^. 

k\ 

2 - For any disjoint subsets A and B, the random variables $a and are 

independent. 

The notion of point process trivially extends to configurations in R'^' x X 
where X is a subset of R™. A configuration is then typically of the form 
{{xn, Un), n > 1} where for each n > 1, a;„ S R'^ and y„ S X. We keep 
writing (a;„, ?/„) as a couple, though it could be thought as an element of R''+™, 
to stress the asymmetry between the spatial coordinate d the so-called 

mark, ?/„. For a marked point process, we denote by <& the set of locations, 
i.e., <i>(ti>) = {X„, n > 1} and by $ the set of both locations and marks, i.e., 
4>(aj) = {{Xn, Yn), n > 1}. A marked point process with position dependent 
marking is a marked point process for which the law of F„ , the mark associated 
to the atom located at Xn, depends only on Xn through a kernel K: 

P(r„ e B I $) = K{Xn, B), for any B d X. 

If X is a probability kernel, i.e., if K{x, X) = 1 for any x G R'^ then it is well 
known that $ is a Poisson process of intensity K{x, dy)dA{x) on R'^ x R™. The 
Campbell formula is a well known and useful formula 

Theorem 1. Let ^ be a marked Poisson process on X R™. Let A be the 
intensity of the underlying Poisson process and K the kernel of the position 
dependent marking. For f : R*"' x R™ ^ R a measurable non-negative function, 
let 

F= / fd^ = Y,f{Xn, Yn). 

Then, 

^[F]= f f(x,y)K{x, dy)dA{x). 

JR'=XR™ 

Definition 2. For F : Fp^fc R, for any x E R'^, we define 

D^F{uj) = F{uj U {x}) - F{lu). 

Note that for F — J /d$, D^F — f{x), for any x S R*^. We now quote from 
0, the main result on which our inequalities are based: 
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Theorem 2 (Concentration inequality). Assume that $ is a Poisson process 
on of intensity A. Let f : R*^ — ^ a measurable non-negative function 
and let 

F{UJ)= J fd<P = J2fiXn{co)). 

n>l 

Assume that \DxF{lj)\ < s for any x £ R'"'. Let 

mF = E[F] = J fix) dAix) 

and 

VF^ J \D,F{Lu)\^dAix)^ J f^ix)dAix). 

Then, for any t G R+, 

^ / „ X f vf f t s 

P(F - mF > i) < exp - g ' 



VF 



where git) ^{l + t) ln(l + t)~t. 



3 Deterministic gain 

We state the following assumptions: 

Assumption 1. The position of each user is independent on the position of 
all other. The users are indistinguishable, i.e., the positions are identically 
distributed. 

Assumption 2. The time between two consecutive demands of users for service 
in the system ( or interarrival time ) is exponentially distributed. 

We define pix) as the surface density of interarrival time in s~^m~^, constant 
in time. Hence, for a region H C B, the mean interarrival rate is h = pix)dx 
in s~^. 

Assumption 3. The service time for every user is exponentially distributed 
with mean 

Assumption 4. The cell C is circular, with radius R and with the antenna in 
the center. 

Assumption 5. The channel gain depends only on the distance from the trans- 
mitting antenna. 

Assumption 6. The surface density of interarrival time is constant. 

These assumptions are commonly done to simplify the mathematical treat- 
ment and are quite reasonable. If we show that the point process given by the 
location of the users is a Poisson process, then it is sufficient to have the two 
first moments in order to apply theorem [2] and then calculate an upper bound 
Psup for Pioss- To do this, we consider the following lemma: 
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Lemma 1. Considering assumptions[^\^and[^ the point process $ of the active 
users positions is, in equilibrium, a Poisson process with intensity dA{x) = 
p{x)v~^dx 

Proof. For a region H , in virtue of assumptions 2 and 3, the number of receiving 
(i.e., active) customers is the same as the number of customers in an M/M/oo 
queue with input rate h and mean service time . It is known Q that the 
distribution of the number of users U in equiUbrium is then 

u'. 

It follows that the first condition of definition[T]is satisfied with intensity measure 



k{H) ^ h/v 



p{x) 



dx. 



H 



Condition 2 of definition [T] follows straightforwardly from assumption [TJ □ 

Without loss of generality, we consider the cell C has its antenna located at 
the origin. We are looking at evaluating 



P{J Nd^> No), 



where N{x) is defined by 



N{x) 



Co 



W\og^ 1 



PtKg 



where g is the mean gain due to shadowing. Note that, with respect to x, N is 
increasing and piecewise constant. Let Rj, j = 1, ■ ■ ■ , N^ax be the values such 
that N{x) — j for x G [Rj, Rj+i)- We can easily determine them by 



PtKg 



According to Theorem [TJ it is then clear that 



1/7 



E 



7Vd4> 



= / A^dA = 



We denote hy rriN the last quantity. Moreover, 
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a 


1.5 


1.6 


1.7 


1.8 


1.9 


2 


p 

^ sup 


0.18 


0.1 


0.04 


0.02 


0.008 


0.003 


A 


0.98 


0.1 


1.15 


1.3 


1.3 


1.4 



Table 1: Comparison between Psup and Pioss for deterministic gain. 



We denote by i^at the last quantity. We take iVo of the form aniN, so that 
according to Theorem [21 



{J N > amN) < Psupia) 



where 

□ / X / VN f {a - l)mNNmc 

'^'"^^"^^^^H'^H 

It is then natural to verify how far this bound is from the exact value of the 
loss probability in simple situations where simulation is available. We used here 
7 = 2.8, Co = 200 kb/s, W = 250 kHz and PtK/I = 1 x 10^. For the surface 
density of interarrival time we use p — 0.0006 min^^m^^ and the service time 
is l/i^ = 1 min, so, the mean number of users in the system is -kR^p/v = 18.85 
users. If we consider the shadowing with a = \/l0 dB and p = 6 dB, we can 
use the mean gain g, giving g = 1/12. Thus, users in the cell boundary use 3 
sub-channels, so Nmax = 3. For a varying from 1 to 2, which corresponds here 
to loss probabilities about 2% or 0.01%, we computed A = \ogiQ Psup/ Pioss- 
Though concentration inequalities are usually thought as almost optimal, the 
results shown in Table [1] seem at first glance disappointing. Remind though 
that the computation of the bound is immediate whereas the simulation on a 
fast PC took several hours to get a decent confidence interval. Remind also that 
the error is about the same order of magnitude as the error made when using a 
usual trick which consists in replacing infinite buffers by finite ones in .Jackson 
networks (see Q). The margin provided by the bounds may be viewed as a 
protection against errors in the modelling or in the estimates of the parameters. 



4 Random gain 



Let us determine now the upper bound probability Psup for Pioss without as- 
sumption[5]but holding all other assumptions of the preceding section. Lemma[l] 
still holds, since it is a consequence of assumptions [H [2] and [3l We also state 
two other natural assumptions: 

Assumption 7. The random gain is totally described by the log-normal shad- 
owing, with mean p and standard deviation a, both in dB. 

For a user at distance d from the origin, the gain is G = 1/5, where S follows 
a log-normal distribution: 



ps{y) 



exp 



2a2 
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where ^ = 10/ In 10. 

Assumption 8. A user is able to receive the signal only if the signal-to-interference 
ratio is above some constant Pmin ■ 

This means, in particular, that the number of subcarriers needed by a trans- 
mitting user is surely bounded by 

The situation is slightly different from that of Section [21 since the functional 
depends on two aleas: positions and gains. Consider now that our configurations 
are of the form (x, s) where x € is still a position and s G R is a gain. Since 
gain and positions are independent, we then have a Poisson process on R'' of 
intensity measure dA{x) ps{y)dy. Thus we want to evaluate an upper bound 
of 

P( J 7Vd$ > No) 



where 



Nix, y) = 



Co 



W\og^ 1 



PtK 
lyxi' 



According to Theorem [21 we must compute 

TTiN = J N{x, y)psiy)dy dA{x) 

and 

VN = sup / \D^^yF{Lj)\'^psiy)dy dA{x) 

N^{x, y)ps{y)dy dA{x). 



Let Po = oo and I3j 
1, let 



1, ■ • • , Nrtiax 



1 , • ■ • , A'^ max 



and Ao — 0. 
Lemma 2. For j = 1, • 

A J = TrR^Q{aj - (Ini?) 



'^{v\\x\\-'<PtK/ipj}Psiy) dy dx 



1. For j 



CxR+ 



7nax ^ f 



+ TTe 



2/C"+2aj/C 



where 



= -(101ogio(Ptif///?, ) - m) and C - ^— r 



Q(Clni?-2/C-a,), 
IO7 
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Proof. We can write 



where Pj = PtK/IjSj. Remind that S is equal in distribution to exp(A/'(/x, tr^)^) 
with ^ = ln(10)/10. Thus after a few manipulations, we get 



where 



Aj = 2Tr I rQ{aj - C In r) dr. 
Q(x) = 



^ / exp(-— )du. 



The final result follows by a tedious but straightforward integration by parts. □ 
Theorem 3. For any function 9 : R — > R, 

9{N{x, y))ps{y)dy dA{x) 

= d{j){Aj-A,-i) + 9{N,na.){nR'-AN„^^^^i). 

Proof. Since N can take only a finite number of values, we have 
9{Nix, y))ps{y)dy dA(a:) 

Y] / ^{(x.y). N{x,y)=j}Psiy)dy dx 

JCxR+ 

Now we see that 

N{x, y)-.7<=^/3j_i<y||xP</3,, 

1 and N{x, y) = N„iax when y||a;p > Pn, 



Nrr 



for j = 1, • • • , N,nax - 

proof is thus complete. 



1. The 
□ 



We used the same set of values as for the simulation of Section [3] together 
with assumptions [8] and [7] with I3min = 0.2. Results of Table [2] show that the 
theoretical bound is rather stable when gains become stochastic. 



a 


1.5 


1.6 


1.7 


1.8 


1.9 


2 


p 

^ sup 


0.2 


0.1 


0.05 


0.02 


0.01 


0.004 


A 


1.7 


1.8 


2.1 


2.3 


2.4 


2.6 



Table 2: Comparison between Pgup and Pioss for random gain. 
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5 General case 



Actually, the method can be applied to more general situations as we illustrate 
now. We consider only assumptions [H El [3] and [H a non-frequency selective 
random gain G with distribution pa, a finite number of antennas with a deter- 
ministic pattern and the assumption that the user will receive the signal from 
the antenna which can provide a better signal-to-interference ratio. 

Now C is the Borel set where it is possible to find users. Possibly, C ~ R^. 
In this region, we have a finite number of antennas J + 1 with the j-th located 
Uj, j = 1, J, and the one we observe located at yo- This means that for each 
user, there is a vector G — (G, G^, G'^) where G is the gain due the antenna 
at yo and due to antenna located at yj. We then define the Poisson point 
process ^' in G x R+~''^, representing the user positions and the gain of each one 
due to each antenna. Again, since gains from different antennas and positions 
are independent altogether, 5' has intensity Xm- 

We define the sets 



Vol 



and 
where 



B = iis,x), s<R{x)), 



Pmin\\x-yoP 

The event {{G,X) £ A') means that the antenna at yo provides the highest 
signal-to-interference ratio to a point X. The event {{S,X) g B) means the 
signal-to-interference ratio provided by the antenna at yo to a point X is higher 
than Pmin ■ By Theorem [21 we are thus led to compute 

iV(||x||,5)"dA„,(x,g), 

AnB 

for = 1, 2. There is no longer a closed form formula for these integrals but 
they can be easily and quickly computed by numerical methods. This yields to 
an upper bound of Pioss ■ We simulated in this section the loss probability for an 
antenna placed at the origin and six other antennas placed at the points yi = 
i2R, 0), ^2 = {R, RV5), 2/3 = i-R, i?V3), yi = (-2i?, 0), = {-R, -RVS) and 
ye — {R, — i?V3), representing an hexagonal arrangement. All other parameters 
are the same as the ones in previous sections. We find a mean tun = 21.60 and 
the second moment vn = 26.81. It turns out that the results are close to the 
results in Section [H suggesting that the approach of Section |4] is satisfactory 
enough with our physical assumptions. 
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6 Concluding remarks 



Using the concentration and deviation inequalities and the difference operator 
on Poisson space, we have calculated the upper bound probability of overloading 
the system by high demand of subcarriers, over path loss and shadow fading. 
To do this we have found the first and second moment of the marked Poisson 
point process of users. We conclude that it is possible to find an upper bound 
for the overloading probability, even in a relatively complex system, which is 
analytically computable in a very simple fashion. The method works for any 
functional of the configurations, possibly enriched by marks, which depends only 
on the positions of each user. It does not work for functionals involving relative 
distance between two or more users. Actually, for such a functional F, there is 
no bound on DxF{uj) valid for all x and uj. 
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